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QCD phase transitions. We discuss the denisty fluctuations which play a significant role to search 
for the critical end point. The importance of spinodal instabilities to distinguish the first-order 
phase transition is also indicated. 
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1. Introduction 



The existence of a critical end point (CEP) in QCD phase diagram is one of the striking expec- 
tations ^ 0], which has been explored based on calculations in effective models and on universality 
arguments. The appearance of the CEP in the temperature T and quark chemical potential pL q plane 
was also studied in terms of Lattice Gauge Theory (LGT) [J3]] . The search for the CEP has re- 
cently attracted considerable attention in the context of heavy ion phenomenology Thus, it is 
of particular interest to identify the position of the CEP in the phase diagram and to study generic 
properties of thermodynamic quantities in its vicinity. The analysis of fluctuations is a powerful 
method for characterizing the thermodynamic properties of a system. Modifications in the magni- 
tude of fluctuations or the corresponding susceptibilities have been suggested as a possible signal 
for deconfmement and chiral symmetry restoration [Q, ||, ||, g). In this context, fluctuations related 
to conserved charges are of particular interest [Qj. 

In this contribution we briefly summarize the properties of the conserved charge fluctuations 
to probe the QCD phase structure. The role of the fluctuations in order to identify the location of 
the CEP as well as the phase boundaries is discussed. We also show that the enhanced baron or 
electric charge density fluctuations could signal the first order phase transition in the presence of 
spinodal decomposition. 



2. Fluctuations of conserved charges 

In our study of fluctuations we adopt the Nambu-Jona-Lasinio (NJL) model as an effective 
chiral model under the mean field approximation [g]. The model describes the chiral phase transi- 
tion where the dynamically generated quark mass M acts as an order parameter. In Fig. |j] we show 
the phase diagram of the two-flavored NJL model for an isosymmetric system in the chiral limit. 
The position of the phase boundary and the tricritical point (TCP) depends crucially on the model 
parameters, like e.g. on the strength of the four-fermion interactions. In the figure we illustrate 
the dependence on the scalar-isoscalar Gs and vector-isoscalar Gy couplings. With increasing Gy, 
the phase transition line at fixed T is shifted to larger \i q due to strong repulsive forces among 
the constituent quarks. Consequently, for sufficiently large value of the vector coupling the TCP 



disappears from the phase diagram [Q, |lfj. 

The position of the phase boundary and the order of the chiral phase transition can be also iden- 
tified through thermodynamic observables, like net baryon number fluctuations which are sensitive 



probes of the phase transition [ ]11| , |12j , |4j, g, g|. Furthermore, fluctuations of conserved charges 
are directly accessible in experiments. Thus, it is of importance to explore the behavior of such 
fluctuations in the vicinity of the phase boundary. 

The quark number and iso-vector susceptibilities, % q and %i respectively, describe the response 
of the net quark density n q and the isovector density n I to the change of the corresponding chemical 
potentials. Thus, % q and %i are defined as derivatives of n q and nj with respect to \i q and [if, 

dn q dnj 

with \i q and \i\ being the net quark and isovector chemical potential respectively. 
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Figure 1: The phase diagram of the NJL model for an isosymmetric system in the chiral limit [foh. The 
tricritical point is indicated by a dot. 



The temperature dependence of % q shows characteristic features, which vary rapidly with \i q . 
The phase boundary is signaled by a discontinuity in the susceptibility. The size of the disconti- 
nuity grows with increasing \L q up to the TCP, where the susceptibility diverges. Beyond the TCP 
the discontinuity is again finite. On the other hand, at \i q = the discontinuity vanishes and the 
susceptibility shows a weaker non-analytic structure at the transition temperature, resulting in a dis- 
continuity of dXq/dT . These critical properties of % q are consistent with that expected for a second 
order phase transition belonging to the universality class of 0(4) spin model in three dimensions 
iilll- 



In Fig. H we show the net quark and isovector susceptibilities along the phase boundary. The 
singularity of % q indicates the existence of the TCP. In the absence of a TCP, the net quark suscepti- 
bility would be a monotonic function of T along the phase boundary, as illustrated in the Fig. |2|-left 
by a dashed-dotted line. We note that on the qualitative level the critical behavior of the net quark 
number susceptibility can be also obtained in the Landau theory |]6[ [HJ: First, the discontinuity 
across the phase boundary vanishes at pL q = 0. Second, the singularity of % q shows up only in the 
chirally broken phase, while the susceptibility in the symmetric phase is monotonous along the 
phase boundary and shows no singular behavior. The isovector fluctuations %h contrary to % q , are 
neither singular nor discontinuous at the chiral phase transition for finite chemical potential. We 
find a rather smooth increase of Xi with increasing jx q along phase boundary line. At the TCP the 
%j remains finite. The non-singular behavior of Xi at the TCP is consistent with the observation 
that there is no mixing between isovector excitations and the isoscalar sigma field due to isospin 
conservation [14]. Recent LGT results [|l^] also show a smooth change of the isovector fluctua- 
tions around deconfinement transition and a fairly weak dependence of Xi on the quark chemical 
potential \i q . 

The net quark number Xq an d the isovector Xi susceptibilities are related with fluctuations of 
the electric charge Xq as 

1 1 1 d 2 P 
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Figure 2: (Left) The net quark number susceptibility % q along the phase boundary in the chiral limit. The 
solid (dashed) line denotes % q in the chirally broken (symmetric) phase. The vertical dotted line indicates 
the position of the TCP. (Right) The isovector susceptibility along the phase boundary [[T~o|] . 
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Figure 3: (Left) The quark number susceptibility near the TCP as a function of the reduced temperature. 
(Right) The critical region where the % q is enhanced by an order of magnitude compared to the free one [|Io|]. 



where P is the thermodynamic pressure. For isospin symmetric system the last term vanishes. 
Hence in this case all relevant susceptibilities are linearly dependent. Clearly, since Xi 1S finite at 
the TCP, the electric charge fluctuations %Q diverge with the same critical behavior as % q . At finite 
\li the properties of %i at th e chiral phase transition change. Then, since the isoscalar sigma field 
mixes with the isospin density [14], the isovector susceptibility exhibits a similar structure as % q , 
with a singularity at the TCP. 

The strength of singularities is governed by the critical exponents whose values are different 
depending on paths approaching the phase transition [15]. The mean-field exponents of the TCP 
and the CEP can be obtained from Landau theory [10]. In Fig. ||]Teft we illustrate the critical 
behavior of % q near the 0(4) critical line and at the TCP. For paths approaching a TCP asymptot- 
ically tangential to the phase boundary, the quark number susceptibility diverges with the critical 
exponent y q = \. On the other hand, approaching the TCP along the first-order transition line, the 
pre-factor of the singular contribution to the quark susceptibility is twice as large as that obtained 
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when approaching TCP along the 0(4) critical line. For other paths the critical exponent is y q = A. 
At the 0(4) critical line, the susceptibility remains finite. The corresponding critical exponent of 
the 0(4) universality class is a ~ —0.2, while in the NJL model under the mean-field approxima- 
tion, a = 0. For non-zero quark mass, at the critical endpoint, the mean-field critical exponent 
along a path not tangential to the phase boundary is 2/3, while along the phase boundary it remains 
equal to unity [@]. When quantum fluctuations are included, the first exponent is renormalized to 
that of the 3D Ising model universality class [16], i.e. e = 0.78. 

In Fig. f|-right we show the "critical" region, where the susceptibility exceeds its value in the 
ideal quark gas by more than an order of magnitude 1 . The differences in values of the critical 
exponents for different "paths" are reflected in the shape of the critical region around TCP. It is 
elongated along the phase boundary, where the singularity is strongest. 



3. Baryon number susceptibility in the presence of spinodal instabilities 

In the previous section we have argued that the enhancement of the baryon number fluctuations 
could be a clear indication for the existence of the critical end point in the QCD phase diagram. 
However, the finite density fluctuations along the first order transition appear under the assumption 
that this transition happen in equilibrium. In non-equilibrium system, a first order phase transition 
is intimately linked with the existence of a convex anomaly in the thermodynamic pressure [|l7|]. 
There is an interval of energy density or baryon number density where the derivative of the pressure, 
dP/dV > 0, is positive. This anomalous behavior characterizes aregion of instability in the (T,n q )- 
plane which is bounded by the spinodal lines, where the pressure derivative with respect to volume 
vanishes. The derivative taken at constant temperature and that taken at constant entropy, 



, and — =0, (3.1) 

dVj T \dVj s 

define the isothermal and isenuopic spinodal lines respectively. 

For finite vale of the quark mass and within large range of parameters the NJL model exhibits 
a critical end point (CEP) that separates the cross over from the first order chiral phase transition. 
The relevant part of the phase diagram in the (T,n q )-plane is shown in Fig. |4| If the first order 
phase transition takes place in equilibrium, there is a coexistence region, which ends at the CEP. 
However, in a non-equilibrium first order phase transition, the system supercools/superheats and, 
if driven sufficiently far from equilibrium, it becomes unstable due to the convex anomaly in the 
thermodynamic pressure. In other words, in the coexistence region there is a range of densities and 
temperatures, bounded by the spinodal lines, where the spatially uniform system is mechanically 
unstable. 

In Fig. |5|4eft we show the evolution of the net quark number fluctuations along a path of fixed 
T = 50 MeV in the (r,n 9 )-plane. When entering the coexistence region, there is a singularity in 
X q that appears when crossing the isothermal spinodal lines and where the fluctuations changes 
the sign. In between the spinodal lines, the susceptibility is negative. Consequently, this implies 



'By "critical" region we mean here the region where the susceptibility is large due to fluctuations and not the region 
where the critical exponents deviates from their mean-field values. 
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Figure 4: The phase diagram of the NJL model jl8[]. The filled point indicates the CEP. The full lines 
starting at the CEP represent boundary of the coexistence region in equilibrium. The dashed curves are the 
isothermal whereas the dotted ones are the isentropic spinodal lines. 



instabilities in the baryon number fluctuations when crossing from meta-stable to unstable mixed 
phase. The above behavior of % q is a direct consequence of the thermodynamics relation 

dP\ nl 1 

w)r-v¥,- <3 - 2) 

Along the isothermal spinodals the pressure derivative in Eq. (3.2) vanishes. Thus, for non- 
vanishing density n q , % q must diverge to satisfy ( |3.2| ). Furthermore, since the pressure derivative 
dP/dV\r changes sign when crossing the spinodal line, there must be a corresponding sign change 
in Xq, as seen in Fig. ||-left. Due to the linear relation between % q , the isovector susceptibility %i an d 



the charge susceptibility Xq fl2.2|), the charge fluctuations are also divergent at the isothermal spin- 
odals. Thus, in heavy-ion collisions, fluctuations of the baryon number and electric charge could 
show enhanced fluctuations across the 1st order transition if the spinodal decomposition appears in 
a system. 

In Fig. ||-right we show the evolution of the singularities from the spinodal lines when ap- 
proaching the CEP. The critical exponent at the isothermal spinodal line is found to be 7= 1/2, 
with Xq ~ (A 4 ~~ l^c) y ', while 7= 2/3 at the CEP [p||. Thus, the singularities at the two spinodal 
lines conspire to yield a somewhat stronger divergence as they join at the CEP. The critical region 



of enhanced susceptibility around the TCP/CEP is fairly small [16, 10], while in the more realistic 
non-equilibrium system one expects fluctuations in a larger region of the phase diagram, i.e., over 
a broader range of beam energies, due to the spinodal instabilities. 

The rate of change in entropy with respect to temperature at constant pressure gives the specific 
heat expressed as 

2s fs\ 2 
Xtt XfiT + — ) Xq ■ 



C P 



1 dT 



TV 



(3.3) 



The entropy Xtt and mixed XfiT susceptibilities exhibit the same behaviors as that of Xq shown in 
Fig. |5|-left. Thus Cp also divergences on the isothermal spinodal lines and becomes negative in the 
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Figure 5: (Left) The net quark number susceptibility at T = 50 MeV as a function of the quark number 
density across the first order phase transition. (Right) The net quark number susceptibility in the stable and 
meta-stable regions [p"8||. 



mixed phase 2 . It was argued that in low energy nuclear collisions the negative specific heat could 
be a signal of the liquid-gas phase transition [19]. Its occurrence has recently been reported as the 
first experimental evidence for such an anomalous behavior 



4. Conclusions 



We presented a brief discussion of probing the QCD phase structure. We have especially 
discussed the importance of conserved charge fluctuations. It was shown that the net baryon number 
susceptibility must yield large contribution around the critical end point (CEP). Consequently, a non 
monotonic behavior of these fluctuations as functions of the collision energy in heavy ion collisions 
could be considered as an indication of the CEP in the QCD phase diagram. 

We have also shown that in the presence of spinodal instabilities the above picture is modified: 
The net quark number fluctuations diverge at the isothermal spinodal lines of the first order chiral 
phase transition. As the system crosses this line, it becomes unstable with respect to spinodal 
decomposition. The unstable region is in principle reachable in non-equilibrium systems that is 
most likely created in heavy ion collisions. Consequently, large fluctuations of baryon and electric 
charge densities are expected not only at the CEP but also when system crosses a non-equilibrium 
first order transition. 
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